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Abstract 

The Old Quantum Mechanics action discretization rules for peri- 
odic motions on the atomic scale (Bohr-Sommerfeld) have been suit- 
ably modified in order to take into account the gravitational field 
instead of the electrostatic one. The new rules are used to calculate a 
few mechanical quantities pertinent to the periodic motions of celes- 
tial objects, and several values are obtained which result in reasonable 
agreement with the corresponding experimental data. 
A gravitational dimensionless structure constant a g is determined, 
using the data relative to the solar system, which allows to quantita- 
tively account for phenomena on a much wider scale. In particular, 
some information is acquired about the recently discovered extrasolar 
planetary systems and about the general empirical law which connects 
the spin of a celestial body with the square of its mass. 
Though no general theory is explicitely proposed in order to explain 
the obtained numerical coincidences, Nelson's stochastic mechanics 
and/or chaotic dynamics seem to be the best candidates for this role. 
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I. Introduction 

In a recent work [?] we suggested that the periodic motions of celestial 
bodies can be described, approximately enough, starting from suitably tuned 
action discretization rules (scaling up of the Bohr-Sommerfeld laws). 
In fact, if the Planck's constant is expressed using the relationship 

e 2 

h = 2,71 , 

a e c 

where a e ~ 1/137 is the dimensionless fine structure constant, the Bohr- 
Sommerfeld rules can be written in the form 



j> Pj dqj = rij 2ir 



e 2 



OLpC 



One can imagine that, if the electron charge was known with higher precision 
in Planck's time, he would have introduced a e instead of h as a new constant. 

The form of the previous equation explicitely shows its 'atomic electro- 
static' purpose. On the other hand, since both the (atomic) electrostatic 
and the gravitational force field follow the r~ 2 law and only differ in their 
'sources': 

e 2 = proton charge times electron charge in the H atom 
GMm = mass m in the gravitational field of a source of mass M , 

we will conjecture that the Old Quantum Mechanics rules hold if one performs 
the substitution 

e 2 GMm 



Qt e (Xg 

where a g is a (new) dimensionless 'gravitational' structure constant to be 
determined. 

Introducing the gravitational potential energy 

GMm 
v , 
r 

it follows by simple calculations that the major semi-axes of the (elliptical) 
orbits are given by 

a n = n 2 ai , 



1 



where n — 1,2,... is the principal number and 

GM 

ai 



a 2 g c 2 



is the 'Bohr radius' of the M gravitational source. 

Other immediate consequences concern the minimum and maximum or- 
bital distances q n \ and Q n \ and the corresponding maximum and minimum 
speeds V n i and v n f. 

f Qn,i = a n (l-e) ; V n j = (l + e)^f 
\ Qn,i = a n (l + e) ; v Hyl = (1 - e)^f ' 

where 



is the elliptical eccentricity, and, as customary, the azimuthal quantum num- 
ber I can take the values 1, . . . , n. The angular momentum is given by 

, GMm 

Jn,l = I • 

OtgC 

In the case of circular orbits (/ = n) one has the constant speed 

a„c 
n 

the corresponding period 



GM „ 
P n = 2ir -3-3 , 



or g c 



and the angular momentum 

GMm 



a g c 



Jn = n . 

Thus, the universal constant 



v * = a g c 



2 



turns out to be the maximum speed for a gravitationally bounded orbitating 
point mass. Correspondingly, 



P* = 2tt 



GM 



which depends on the strength M of the gravitational source, turns out to 
be its minimum period, and 



GMm 



OtgC 



its minimum angular momentum. 



II. A 'domestic' test case: orbital motions around the Sun 

As a first test case we can apply the previous formulae to the motion of 
the solar planets. In this case we use the following relationships 



Major semi-axes : a n 



n 2 af 



where a{ 







Periods : 



Mean speeds : 



P n = n 3 P? , where Pf 







GM® 



2nGM e 



where M® indicates the mass of the Sun. 

In order to determine a g , we previously need to select the quantum prin- 
cipal numbers for the various planets. Note that, when a given number for 
a given planet is chosen, the numbers for all the other planets turn out to 
be determined by their observed mean speeds. Of course, the results are 
significant if the concerned quantum principal numbers are small, and if no 
satisfactory agreement exists for semi-integer n values. 

One can easily verify that the only consistent choice lies in assigning to 
Mercury the principal number 3. Using the corresponding known speed one 
obtains the initial guess for the gravitational structure constant 



ot g = 3 



^Mercury 



and the following quantum number selection: Mercury^3, Venus— >4, Earth^5, 
Mars— >6, Jupiter— Saturn^l5, Uranus— >-21, Neptune— >26, Pluto— >30. 
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On this basis, one can use the relationship v n = a g c/n to refine the 
estimate of a g through a general least square fit. 
In this way, the following results can be obtained: 

1/ag = 2086 ± 14 , 

= 143.7 ± 1.0 km/s , 
af = 0.04297 a.u. , 
Ff = 3.269 d , 

and 



Planet 


n 


v bs (km/s) 


Vcaic (km/s) 


Mercury 


3 


47.87 


47.90 


Venus 


4 


35.02 


35.92 


Earth 


5 


29.78 


28.74 


Mars 


6 


24.13 


23.95 


Jupiter 


11 


13.06 


13.06 


Saturn 


15 


9.64 


9.58 


Uranus 


21 


6.80 


6.84 


Neptune 


26 


5.43 


5.53 


Pluto 


30 


4.74 


4.79 



In Fig. 1 the v n = a g c/n and a n = n 2 a\ relationships are graphically 
compared with the actual orbital mean speeds and the actual major semi- 
axes of the planets of the Sun. One can see that the solar planetary motions 
are reasonably accounted by our model of 'gravitational atom', even for small 
principal numbers. 

As the azimuthal numbers are concerned, considering one of the two plan- 
ets with valuable eccentricity, i.e. Pluto (n = 30), and using the eccentricity 
formula with / = 29, one obtains e ca i c = 0.256, in fair agreement with the 
observed value e b s ~ 0.25. Unfortunately, no allowable calculated value fits 
satisfactorily the eccentricity of Mercury (whose motion shows on the other 
hand the well known anomalies). 

Very recently [?] a new object has been discovered, named 1996TL 66 , 
which can be assumed to adequately represent a number of bodies located 
between the Kuiper belt (beyond the Neptune orbit) and the Oort cloud 
(which reachs abouto 50,000 a.u. from the Sun). The orbit of this object, 
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the brightest trans-neptunian body but Pluto and Charon, shows the very 
large eccentricity e = 0.58 and a major semiaxis a n = 83.7a.-u.. 

Using our mechanical model, one gets the (n,a n ) pairs (43,79.4 a.u), 
(44,83.2 a.u.) and (45,87.0 a.u.): thus, the principal number n = 44 seems 
to be adequate to this 'new planet'. 'Scanning now the allowable principal- 
azimuthal number pairs we obtain (44,35) — > e = 0.61, (44,36) — > e = 0.58 
and (44, 37) — > e = 0.54, so that the pair (n, /) = (44, 36) can be assigned to 
1996TL 66 . Nevertheless, we explicitely agree on the fact that results based 
on such large quantum numbers can not contribute to confirm or disprove 
our model. 

III. Orbital motions around the planets 

On the basis of the assumed universality of the constant a g one can use, in 
evaluating the 'Bohr radius' of each planet, the formula 

planet) _ 

Through r ^ lanet ^ one can easily calculate the allowed orbits relative to each 
planet. In particular, one can determine the principal number rif which cor- 
responds to the first actual 'free' orbit (i.e., to the allowed orbit immediately 
out of the planet body). In the following table we show some interesting 
results (all the lengths are reported in km, and the general relative uncer- 
tainty of the calculated radii is about 1%). In the last two columns we report 
the greatest 'not-free' calculated orbital radius and the observed equatorial 
radius of the concerned body. Note that, in estimating n/ ree , the well known 
internal stability Roche limit [?] has not been taken into account: some un- 
derestimation could result. 

Considering the principal number of the first 'free' orbit, one can note 
that that for the internal planets (and for Pluto) these orbits correspond to 
relatively high (> 19) principal numbers. The opposite, i.e. small (< 10) 
'free' quantum nubers, is verified for the big external planets. 
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Bodv 


Bohr radius 


77, /■ „ 
' v jree 


'v free 


J* , 

(L free 1 


Eq. radius 


Sun 


6.50 10 6 


l 


(6.50 ±0.06) 10 6 




6.96 10 5 


Mercury 


1.08 


48 


2488 


2386 


2439 


Venus 


15.86 


20 


6344 


5725 


6051 


Earth 


19.45 


19 


7021 


6302 


6378 


Mars 


2.09 


41 


3513 


3344 


3397 


Jupiter 


6190 


4 


99040 


55710 


71492 


Saturn 


1853 


6 


66708 


46325 


60268 


Uranus 


282 


10 


28200 


22842 


25559 


Neptune 


333 


9 


26918 


21312 


24764 


Pluto 


4400 10~ 4 


160 


1126 


1112 


1123 



If the orbits of the natural satellites of the planets are considered one 
really obtains ambiguous results. For instance, only three of the four greater 
Jupiter satellites confirm our rule. This fact might denote some failure of our 
model or, on the contrary, that the fourth satellite does not originate from 
the same Jupiter-material as its companions. A particular case is that of the 
Earth-Moon pair: we will consider this case later. 

As the Saturn rings are concerned, one obtains the following results, where 
the unit lenght is 10 3 A;m, and IE, OE indicate respectively the internal and 
the external edge of the considered ring: 



Rings 


D IE 


ClE 


BlE 


Cassini 


A IE 


F 


G 


ElE 


Eoe 


^obs 


67.0 


74.5 


92.0 


119.8 


122.2 


140.4 


170 


180 


480 


n 


6 




7 




8 


9 




10 


16 


T'calc 


66.7 




90.8 




118.6 


150.1 




185.3 


474.4 



The uncertainties of the calculated values are all about 1%. The calculated 
inner edge positions of most rings are in fair agreement with the correspond- 
ing experimental values. Note that the boundary of the C ring is presum- 
ably determined by the interactions between the charged dust particles and 
the planetary magnetic field, so that a purely gravitational theory does not 
account for this fact. Note also the remarkable coincidence between the cal- 
culated and the observed values of the boundaries of the broad rarefied E 
ring. 
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IV. Extrasolar planets 

During the last two years a few extrasolar planetary systems have been 
discovered on the basis of the detected periodical anomalies in the motions 
of their stars. Unfortunately, the mass of most the concerned stars are at 
present estimated with some inaccuracy. Almost all the discovered planets, 
with masses ranging between one half and ten times the Jupiter mass, turn 
out to move somewhat near to the star, their guessed major semi-axes ranging 
from about 0.05 a.u. and 1 a.u. ( only three planets show a mean distance 
from the central star greater than 1 a.u. and less than 2 a.u.). 

In order to test our scheme, we have considered only the 19 extra-solar 
planets which have been confirmed (at present) [?] and two more planetary 
systems (Proxima Centauri and Barnard's) whose star masses are known 
with sufficient accuracy. The star masses of the other 19 confirmed plan- 
etary systems have been guessed from the star type information using the 
values M st reported in literature ([?]). Then, the major semi-axes a st of the 
planetary orbits have been calculated from the corresponding measured pe- 
riods. Unfortunately, the starting point (i.e., the values of the star masses) 
was somewhat inaccurate. 

Thus, an alternative path has been devised: starting from the measured 
planetary periods, and using the relationship 

we scanned the various hypothetical star masses Mi, M 2 , M n which 
would correspond to the various planet principal numbers 1,2, ... ,n. Then, 
the mass value Mj nearest to the guessed star-type mass M st was singled 
out: thus, the the appropriate principal number j was determined , and this 
in turn allowed an appropriate estimate a^ st ^ of the major semi-axis. In the 
following table, in order to put in evidence the ambiguity level in choosing 
the 'most likely' star mass, we report, together the chosen mass Mj, its 
neighbouring hypothetical values Mj_i and Mj +1 . 

We explicitely note that for most of the considered planetary systems the 
choice of the appropriate Mj is unequivocal enough. 

For the last two planetary systems, the planet period was astrometrically 
determined. In both cases, the chosen Mj value turns out to fit exactly the 
measured one: this fact seems to support the correctness of our estimation 
method. 
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Star 


Type 


M st 


p(obs) 




M,-_i 


M 7 - 

J 


M,-+i 


n 


-(est) 
a„ 


51 Peg 


G2IVA 


1.05 


4.229 


0.050 




1.30 


0.16 


1 


0.056 


ups Androm. 


F7V 


1.20 


4.611 


0.057 




1.42 


0.18 


1 


0.061 


55 Cancer 


G8V 


0.90 


14.648 


0.110 


4.50 


0.56 


0.17 


2 


0.097 


rho CrB 


G0V-G2V 


1.05 


39.645 


0.230 


12.18 


1.52 


0.45 


2 


0.261 


16 Cyg B 

J o 


G2.5V 


1.05 


804.000 


1.720 


1.98 


1.14 


0.72 


6 


1.766 


47 Uma 


GOV 


1.05 


1088.445 


2.110 


1.55 


0.98 


0.65 


7 


2.050 


tau Bootis 


F6IV 


1.30 


3.313 


0.046 




1.02 


0.13 


1 


0.044 


70 Virgo 


G4V 


0.90 


116.600 


0.430 


4.48 


1.33 


0.56 


3 


0.512 


HD 114762 


F9V 


1.20 


84.050 


0.300 


3.23 


0.96 


0.40 


3 


0.369 


HD 110833 


K3V 


0.73 


270.040 


0.800 


1.30 


0.66 


0.38 


5 


0.712 


BD-04 782 


K5V 


0.67 


240.920 


0.700 


1.16 


0.59 


0.34 


5 


0.635 


HD 112758 


K0V 


0.79 


103.220 


0.350 


1.17 


0.50 


0.25 


4 


0.340 


HD 98230 


F8.5V 


1.30 


3.980 


0.060 




1.22 


0.15 


1 


0.052 


HD 18445 


K2V 


0.73 


554.670 


0.900 


1.36 


0.79 


0.50 


6 


1.219 


HD 29587 


G2V 


0.98 


1157.843 


2.500 


1.65 


1.04 


0.69 


7 


2.180 


HD 140913 


GOV 


1.05 


147.940 


0.540 


1.68 


0.71 


0.36 


4 


0.488 


HD 283750 


K2 


0.73 


1.790 


0.040 




0.55 


0.07 


1 


0.024 


HD 217580 


K4V 


0.70 


454.660 


1.000 


1.12 


0.65 


0.41 


6 


0.999 


Alpha Tau 


K5III 


1.20 


654.000 


1.350 


1.61 


0.93 


0.59 


6 


1.437 


Prox.Ccnt. 




0.10 


42 




0.20 


0.10 


0.06 


5 


0.111 


Barnard's 




0.12 


132.0 




0.19 


0.12 




7 


0.249 



V. Quantized galactic redshifts 

Starting from the seventies, various authors ([?], [?], [?],[?], [?]) found 
that the frequency analysis of the galactic redshifts (both for binary galaxies 
and in the general case) indicates speed differences which are multiples of 
37.5 km/s (Tifft suggests 36 km/s). Using our model, one can try to explain 
both this discretization and the estimated speed periodical structure as effects 
of the speed quantization rule 

a„c 143.7 , . 

v n = - JL - = km/s . 

n n 



8 



VI. Spin of the celestial bodies 

In order to validate or disprove the effectiveness of our scheme and the 
universality of the calculated a g , we have taken consideration the known 
empirical the relationship between the angular momentum J of a celestial 
body and its mass M. Studying this question, both from a theoretical and 
a phenomenological point of view, Brosche (see for instance [?], [?], [?]) and 
later Wesson ([?], [?]) suggested that the law 

J = pM 2 

holds for purely gravitationally bounded systems. Wesson, fitting the empir- 
ical values and taking into account some theoretical consistency constraints 
suggested to use the value 

p ~ 8 KT 16 cm 2 g- 1 s~\ 

Consistently with the correspondence e 2 < — > GMm on which our me- 
chanical model is based, and recalling that the electron spin can be written 
in the form 

2 a e c ' 

we can suppose that the rotational momentum of a purely gravitationally 
bound celestial body is given by 

_ 1 GM 2 

2 a g c 

As a consequence, we get p = G/(2a g c) = 2.32 1CT 15 cm 2 g^ 1 s^ 1 . 

We note that, independently from this new guess, the previous formula 
can also be directly obtained from our considerations on the discretization of 
the orbital motions if a flat continuous purely gravitating disk is considered in 
which each infinitesimal mass ring rotates on the Bohr radius of the internal 
part of the disk which it delimits. In this case, the total angular momentum 
can be easily calculated by the formula 



/G r 
dJ = / m dm , 
a g c J 



which just gives the previously guessed expression. In fact, this remark sug- 
gests that the concerned expression works well for flat or almost flat celestial 
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bodies (prevailing gravitational energy) whereas some correction factor can 
be needed in the opposite cases (when both the gravitational and the elec- 
trical energy play their role). 

Anyway, we have compared these theoretical spin values with empirical 
values found in the literature or, as in the case of most planets, with the 
values 

2 2 2vr 

Jobs — ^MR — — , 
obs 

where P b s indicates the observed rotation period (note that this formula 
refers to spherical and homogeneous rigid bodies and in fact it gives upper 
limits for the actual values). 

In Fig. 2 we report the logarithmic graph of our expression J = pM 2 com- 
pared with the experimental data: since no fit was performed, the agreement 
between our formula and the data (over more than twenty orders of mass 
magnitude) can be judged quite good. 

In fact, the comparison shows some small discrepancy for the the so- 
lar planets, which on the other hand are examples of celestial bodies held 
together by other forces besides the gravitational one. 

In the following table the calculated and the "observed" (upper limits) 
angular momenta for some planets are shown: the Earth apart, only for 
Saturn, owing to its oblate shape, the angular momentum J = 8.16 1 37 J s 
can be directly calculated without using the previously remarked assumptions 
about the spherical planet form and its mass distribution. 
A clear discrepancy regards the Earth, which shows J ca i c = 8.24 10 33 Js to 
be compared with = 5.88 10 33 J s. This fact can probably be accounted 
by the permanent loss of angular momentum caused by the tidal friction. An 
increasing rate of the rotation period of the Earth of about 16 s/10 6 years ~ 
3 10~ 12 has been estimated [?]. This would correspond to a Moon braking 
action beginning 'only' about 1600 thousand millions of years ago, long after 
the Earth ocean was formed. On the other hand, the discrepancy in question 
between calculated and measured angular momentum of the Earth would 
be even greater if the whole Earth-Moon system ( J obs = 3.47 10 34 J s) was 
considered. Both this and the previous remarks seem to indicate that the 
Moon and the Earth can not be considered as a single celestial body: both 
could have originated from protoplanetary material on the same orbital ring 
(n = 5, as ~ 1.07 a.u.) and subsequently have gravitationally captured each 
the other (actually the Earth captured the Moon) about 2 10 9 years ago. 
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This hypothesis would also account for the 9% of discrepancy between the 
observed and calculated terrestrial orbital radius round the Sun. 



Body 


Jcalc 


[Js] 


J obs 


[Js] 


Mercury 


2.53 


10 31 


< 1.02 


1Q 30 


Venus 


5.50 


10 33 


< 2.14 


10 31 


Earth 


8.24 


10 33 


5.88 


10 33 


Mars 


9.56 


10 31 


< 2.10 


10 32 


Jupiter 


8.37 


1Q 38 


< 6.83 


10 38 


Saturn 


7.50 


10 37 


8.10 


10 37 


Uranus 


1.76 


10 36 


< 2.50 


10 36 


Neptune 


2.46 


10 36 


< 2.30 


10 36 


Pluto 


5.23 


10 28 


< 1.5 


10 29 



VII. Conclusions 

• Assuming that the previously stated analogies hold between celestial 
and Old Quantum Mechanichs, one can calculate values for a number 
of celestial quantities which turn out to be in remarkable agreement 
with the corresponding observed values. 

• The gravitational structure constant a g ~ 1/2086 has been used in 
order to account some observed values related to 

— the planetary motions around the Sun 

— the structure of the Saturn rings 

— the extra-solar planetary systems 

— the quantization of the galactic red-shifts 

— the spin of a broad class of celestial bodies, from solar planets to 
supercluster. 

• No definite theory has yet been devised on the why the concerned 
analogies hold. Nevertheless, Nelson's stochastic mechanics ([?], [?]) 
seems at present the best candidate to offer valuable suggestions. In 
fact, in Nelson's scheme, a particle of mass which moves in a potential 
field is always subjected to a Brownian motion with diffusion coefficient 
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h/2m: thus, the particle motion can be described using a probability 
density which turns out to obey the Schrodinger equation (of which 
the rules of the Old Quantum mechanics are consequences). If one 
consider the gravitational potential of a source of mass M, the diffusion 
coefficient of the appropriate Brownian motion would be GM/(2a g c). 
This approach has been already applied to the study of the solar system, 
with some interesting results ([?], [?]). 

In fact, we suspect that the Nelson's stochastic approach is an effective 
way to describe a purely classical deterministic chaotic dynamics, where 
the chaotic behaviour arises from the strong nonlinearities in a many 
body problem. The verified existence of stable attracting orbits in 
several many body problems can be a clue to the truth of this hypotesis. 

A different interesting approach, which gives results very similar to 
those previously mentioned, has been worked out by Nottale (see for 
instance [?],[?]). It is based on a scale relativistic scheme, where the 
space-time resolution is an essential physical variable, and the physical 
laws are to be scale-covariants under resolution transformations: the 
fractal structure of the space-time plays a role of capital importance in 
the theory. 
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